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The IceCube detector took data in its 22-string configuration in 2007-2008. This data has been analyzed to 
search for extraterrestrial point sources of neutrinos using several methods. Two main methods are discussed and 
compared here: the binned and the unbinned maximum likelihood method. The best sky-averaged sensitivity 
(90% C.L.) is E^^,,^ = 1.3 X 10"" TeVcm"^s"^ to a generic E"^ fiux of i^^ over the energy range from 3 TeV 
to 3 PeV. No neutrino point sources are found from the individual directions of a pre-selected catalogue nor in 
a search extended to the northern sky. Limits are improved by a factor of two compared to the total statistics 
collected with the AMANDA-II detector and represent the best results to date. 



1. Introduction 

IceCube [ [T] is a neutrino telescope located at 
the South Pole. Its major goal is to discover 
high energy neutrinos of extraterrestrial origin. 
Its vertical strings, each equiped with 60 optical 
modules (OMs), are being deployed year by year 
in the ice at depths between 1.45 and 2.45 km. 
In the 2007/08 season 22 strings (IC-22) were 
operating. Results from this dataset are given 
here. By 2011 IceCube will be completed with 80 
strings extending to 1 km'' instrumented volume. 

2. Methods and Event Selection 

Two independent analyses developed for point 
source searches in the TeV - PeV energy region 
are presented here: a binned method, an event 
counting approach using a fixed circular angu- 
lar search bin; and an unbinned method, a likeli- 
hood approach based on the point spread function 
and a simple energy estimator. In a declination 
band, a negligible point source neutrino contri- 
bution compared to the total background is as- 
sumed, allowing an estimate of the background 
from scrambling data. 

Muon tracks are reconstructed using likeli- 
hood algorithms, taking into account recorded 
hit times and charges. Neutrino event selection 



for both analyses is based on similar parameters: 
angular cuts, which partially filter out southern 
sky atmospheric muons {^atm), and track qual- 
ity parameters, which help further reject mis- 
reconstructed fJ-atm- These quality parameters in- 
clude angular uncertainty estimates, the unscat- 
tered photon hits with a small time residual w.r.t. 
the Cherenkov cone, a likelihood ratio comparing 
a prior of a fiatm track to the reconstructed one, 
a reduced likelihood of the directional reconstruc- 
tion, etc. The main source of residual background 
is atmospheric neutrinos v^i^i^- 

2.1. Binned search 

The binned method [ [2] distinguishes a local- 
ized excess of signal from a uniform background 
(^fj-atm ^^'^ misreconstructed Hatm) using a circu- 
lar angular search bin. The sensitivity, in a given 
declination band, is calculated by comparing the 
number of simulated signal events reconstructed 
inside the search bin, risig , to the average number 
of background events from data, using: 



(ribg) 



(1) 
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*See |http://www. icecube. wise. edu| for a full list of authors. 



where $0 is the signal fiux normalization, FC^g^ , 
given by Poisson statistics in the Feldman & 
Cousins approach ['3^ represents the average event 
upper limit (90% C.L.) for the expected back- 
ground in the bin, Ubg, and no true signal. 
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Figure 1. Left: IC-22 neutrino effective area in declination ranges (binned method). Right: IC-22 
sensitivity and discovery potential to point sources with E~'^ spectrum as functions of declination. 



Cuts on track quality parameters are optimized 
considering the variation in detection and recon- 
struction efficiencies with declination. In addi- 
tion, the search bin radius depends on declination 
(mean value = 2.1°). Final cuts were chosen so as 
to optimize sensitivity for a compromise between 
E^'^ and neutrino energy spectra. 

2.2. Unbinned maximum likelihood search 

The unbinned maximum likelihood method [ 
|4] constructs a likelihood function, C, which de- 
pends on the signal probability density function 
(pdf) and background pdf, for a given source lo- 
cation, Zs, and total number of data events ritot- 

"tot / \ 

C=T\ (^S^ix^,Xs,E,,^) + ^B,ix^,E,)) .(2) 
r^-^ \ntot ntot J 

The signal pdf. Si, is characterized by each 
event's reconstructed direction, gular uncer- 
tainty around Xs, energy estimation, Ei, i.e. the 
number of OMs hit, and 7, the assumed energy 
spectrum index {E~'^). The background pdf, Bi, 
is determined from the declination distribution 
and energy estimators of the data sample. The 
number of signal events, ris, is found by maximiz- 
ing the likelihood ratio of the background plus sig- 
nal hypothesis against the background-only case. 

This method does not need tight cuts, as in the 
binned analysis, because the likelihood incorpo- 
rates energy information that can separate signal 
from background. A harder event selection does 
not further improve the final sensitivity. 



3. Data Samples 

The declination range from -5° to 85° has been 
scanned in the TeV-PeV energy region using IC- 
22 data with a lifetime of 276 days. Each analysis 
has followed its own event selection criteria, ar- 
riving at a final sample of 5114 (2956) events for 
the unbinned (binned) method. Since the binned 
method needs more efficient cuts to get a bet- 
ter sensitivity, the final sample is smaller than 
the unbinned method, which contains more low 
energy t'/j^j^. The content of the sam- 

ples is greater than 90%, where the contamination 
comes from misreconstructed down-going muons, 
mostly time-coincident multiple events. 

From simulation, a sky-averaged median angu- 
lar resolution of 1.4° is estimated for signal neu- 
trinos with spectrum. The final effective area 
for a lyfj^+Pfi flux is presented in Fig. [l]and is simi- 
lar for both analyses. The effect of the absorption 
in the earth, due to increased cross section with 
energy, can be seen for events at large declination. 

4. Results 

Sensitivities and discovery potentially, as func- 
tions of declination, are compared in Fig. [T] The 
sky-average senstitivity at 90% C.L. to a generic 
E-2 flux of is 1.3 (2.0) xlO~" TeV^icm-^s-i 
(E/TeV)^^ for the unbinned (binned) method. 
The enhancement, without including systcm- 

-"^ Smallest flux detectable at 5ct in 50% of trials. 
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atic uncertainties, gained by using the unbinned 
method is 35%. The unbinned sensitivity shows 
a factor of 2 improvement over that of the total 
statistics collected by AMANDA-II [ [S] and rep- 
resent the best limits to date. 

A pre-defined list of 28 sources (galactic and ex- 
tragalactic) gave no significant excess over back- 
ground. The major deviation was found for lES 
1959-1-650 (Crab Nebula) with a post-trial p- value 
of 66% (91%) by the unbinned (binned) analysis. 

The results from the all-sky search for both 
analyses are shown in Fig. [21 The maximum 
upward deviation from the background hypoth- 
esis is located at r.a. — 153° (183.5°) , dec 
= 11° (47°), with a post-trial p-value of 0.67% 
(67.6%) for the unbinned (binned) analysis. The 
final corrected p-value for the unbinned analysis 
is 1.3%. This represents the overall outcome of 
IC-22, as it was decided a priori that the final 
result would be based on the more sensitive un- 
binned analysis. The unbinned-spot, to be ver- 
ified with IC-40 data, was further investigated 
(e.g. time-dependent analysis) resulting in no 
significant event cluster in time (no contributing 
events are closer together than 10 days). 



5. Summary & Outlook 

The results of both methods are consistent 
with background fluctuations. No neutrino point 
sources are found from the individual directions 
of the pre-selected catalogue nor in the extended 
northern sky search. The limits set with IC-22 are 
the best to date. Systematic uncertainties will be 
included in the final results. 

IceCube is half completed with 40 strings tak- 
ing data since April 08. A better sensitivity is ex- 
pected given the bigger effective area and better 
angular resolution. In the coming 2008-09 season, 
16 more strings are planned to be deployed. 
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